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ABSTRACT 

Based on MAGIC observations from June and July 2007, we have obtained an integral upper limit to 
the VHE energy emission of the globular cluster M13 of F(E > 200 GeV) < 5.1 x 10~ 12 cm~ 2 s~\ 
and differential upper limits for E > 140 GeV. Those limits allow us to constrain the population of 
millisecond pulsars within M 1 3 and to test models for acceleration of leptons inside their magnetospheres 
and surrounding. We conclude that in Ml 3 either millisecond pulsars are fewer than expected or they 
accelerate leptons less efficiently than predicted. 

Subject headings: gamma rays: observations, globular cluster: individual (M13) 
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1. Introduction 

Globular clusters (GC) are very interesting sites for 
probing high energy processes due to their large con- 
tent of evolved objects. Millisecond pulsars (MSP) 
constitute a large fraction of these objects and it has 
been estimated that a ty pical massive GC contains of 
the order of 100 of them (lTavanilll993l) . Moreover, the 
largest sample of MSPs discovered up to now in radio 
observations are located in th e GC Ter 5 (23 MSP) , 
and Tuc 47 (22 MSP) (see e.g. lCamilo & Rasioll2005k . 

Fluxes of TeV 7-rays from GC detectable by cur- 
rent Cherenkov telescopes have been predicted based 
on estimates on the population of MSPs and the effi- 
ciency of lepton acc eleration in th eir surrounding (se e 
Bednarek & Sitarekl d2007l) and IVenter el al.l d2009l) ). 
These 7-rays would be produced by accelerated lep- 
tons scattering off photons of the microwave back- 
ground radiation or the thermal emission of an ex- 
tremely dense cluster of solar mass stars inside the 
GC. Acceleration of leptons could take place in i) the 
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shocks within the GC, coming from the collision of 
the winds of MSPs, or ii) in the pulsar inner mag- 
netosphere or their wind regions. In addition, 7- 
rays in the sub-TeV energy range could also be orig- 
inated in the inner MSP magnetos phere directly, as it 
is predicted in the c alculations by Buli k et al.l (2000); 
lHarding. Usov & Muslinovl d2005l) . or could be pro- 
duced in the vicinity of radio emitting bloc ked pulsars 
dAharonian et alJl2005[ lAlbert et alJl200fjh inside low 
mass binary systems dTavanilll99ll) . 

GCs have been observed occasionally by Cherenkov 
Telescopes to probe for this possible VHE 7-ray emis- 
sion. The few experimental results reported in the 
literature are upper limits on t he emission of M l 3 by 
the WHIPPLE Collaboration dHall et al 1 120031). M1 5 
by the VERITAS Collaboration dLeBohec et al.ll2003l) . 
and u> Centauri by the CANGAROO Collaboration 
dKabuki et alJ 12007). Very recently, the Fermi LAT 
telescope has detected high energy 7-ray emission 
(E > 100 MeV) from one of the closes t and most mas- 
sive GC, Tuc 47 dGuillemot et alj2009l) . and HESS has 
obtained an upper limit of 6.7 x 10~ 13 ph cm" 2 s 



for energies E > 800 GeV dAharonian et al J 12009) 



but given the possible complexity of the emission 
in the GeV range, it is not possible to establish any 
connection between these results. This HESS result 
constrains the magnetic field in the pulsar nebula as 
a functio n of the number of MSP in the GC for the 
model by IVenter el al. (2009), and in the efficiency of 
the rotational energy conversion of MSPs into rela- 
tivistic leptons for the model by Bed narek & Sitarekl 
d2007l) . 

In this paper we report the results of observations 
with the MAGIC telescope of the globular cluster 
M13, and we discuss the constraints that our results 
impose to the population of millisecond pulsars and 
their lepton acceleration efficiency. Ml 3 belongs to 
the class of normal globular clusters, and has an esti- 
mated mass of 6 x 10 5 M Q . It is located in the North- 
ern sky at a distance of 7 kpc. Its core radius i s about 
~ 1.6 pc, with a half mass radius of ~ 3.05 pc ( Harris! 
19991) . By now 5 millisecond pulsars have been de- 
tected in Ml 3, with periods ranging between 2 and 10 
ms. The aforementioned observation of this cluster in 
search f or VHE emissio n by the WHIPPLE Collab- 
oration dHall et al.l 12003) led to a flux upper limit of 
1.08 x 10" 11 ph. cm" 2 s" 1 at energies E > 500 GeV. 
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2. Observations and data analysis 

The MAGIC telescope is an Imaging Atmospheric 
Cherenkov Telescope (IACT) located at the Observa- 
tory Roque de los Muchachos on the Canarian Island 
La Palma. It has an exceptional light detection ef- 
ficiency provided by the combination of a 17 m di- 
ameter mirror and a pixelized camera composed of 
576 high quantum efficiency, hemispherical photomul- 
tiplier tubes (PMT). This allows MAGIC to reach a 
standard trigger threshold of ~ 60 GeV. For energies 
above 150 GeV, angular and energy resolutions of the 
telescope are ~ 0.1 deg. and ~ 25% respectively (see 
lAlbertetalJ E008a for further details). Besides this, 
in February 2007 its data acquisition system was up- 
graded with multiplexed 2 GHz Flash Analog-Digital 
converters which improved the timing resolution of the 
recorded shower images. Accordi ngly the sensitiv ity 
of MAGIC improved significantly dAliu et al.ll2009h to 
1.6% of the Crab Nebula flux above 270 GeV for 50 
hours of observation. 

We observed Ml 3 at zenith angles ranging from 
8° to 31° between June 12th and July 1 8th of 2007 



in fal se-source tracking (wobble) mode (iFomin et al. 



11994 . with two directions at 24' distance and oppo- 
site sides of the source direction. This technique al- 
lows for a reliable estimation of the background with 
no need of extra observation time. The collected data 
amount to 20.7 hours after rejecting events affected by 
unstable hardware or environmental conditions. Be- 
sides this, events with a collected charge below 300 
photo-electrons were rejected in order to maximize the 
analysis sensitivity. This selection resulted in a sample 
with a peak energy of 190 GeV. 

Data analysis was carried out using the standard 
MAGIC analysis and reconstruction software chain, 
which proceeds in several steps. Initially, a standard 
calibration of the P MT signal pulses is performed 
^Albert et al.1 l2008bl) . Then, pixels containing no 
useful information for the shower image reconstruc- 
tion are discarde d by an image cleaning procedure 
dAliu et al.l 120091) . Afterwards e vent image parame- 
ters are calculated ( Hil lasl 1 1 985D using the surviving 
pixels. In addition to the classical Hillas parameters, 
two timing parameters are computed, namely: the gra- 
dient of the arrival times of the Cherenkov photons 
along the shower axis; and their arrival time spread 
over the whole shower. The signal-to-noise maxi- 
mization is achieved using a multidimensional clas- 
sification procedure based on the Random Forest (RF) 



method dAlbert et alj|2008cl) . where a hadron likeness 
measure, the so called hadronness, is computed for 
each event based on the image and time parameters. 
Moreover, a regression RF trained with a Monte Carlo 
simulated 7-ray sample is used to estimate the energy 
on an event by event basis. Finally the angle between 
the major axis of the shower image ellipse and the 
source position in the camera, the so called Alpha an- 
gle, is used to select 7-ray candidates in the direction 
of the source. To estimate the residual background, the 
angle Alpha is also computed with respect to a posi- 
tion symmetric to the source position with respect to 
the camera center. In what follows, this position used 
to estimate the residual background is referred as the 
background region. 

Main contributions to the systematic error of our 
analysis are the uncertainties in the atmospheric trans- 
mission, the reflectivity of the mirros (including losses 
due to surface roughness) and the light catchers, the 
photon to photo-electron conversion calibration and 
the photo-electron collection efficiency in the photo- 
multiplier front-end. A detailed discussion of their 
contribution to the flux uncertainties can be found in 
lAlbert et alJ ((2008a), where they are estimated to add 
up to 30% of the measured flux value. 

3. Results 

Figure Q] shows the obtained distribution of the Al- 
pha angle for the source region and the estimated back- 
ground. It has been obtained for events surviving a 
hadronness cut tuned to yield an energy independent 
7-ray selection efficiency of 80%, estimated by means 
of a Monte Carlo simulation. We define the signal re- 
gion as the smaller interval in Alpha angle that con- 
tains the 80% of the 7-rays for each energy bin, es- 
timated using a Monte Carlo simulation. Their lower 
bounds are at Alpha= and the upper ones are shown 
in second column of Table Q] for each energy bin. We 
find a total of — 23 ± 57 excess events after background 
subtraction in this signal region for an energy range 
extending from 140 GeV to 1.1 TeV. In addition a 
search for signals in a region of 1 deg of radius around 
Ml 3 yields no positive detection. We have obtained 
upper limits to the VHE flux from Ml 3 for different 
energy bins, shown in Table [T] Th ese have been com- 
puted using the Rolke method by IRolke et al.l {2005) 
at a 95% confidence level, and they take into account 
a 30% of systematic uncertainties in the flux level. 
The upper limit to the integral flux for energies above 
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E = 200 GeV, assuming a spectral index of 2.6, is 

5.1 x 10~ 12 cm- 2 s" 1 . 

Table 1 : Differential upper limits 



Energy bin 


Upper Alpha 


Events 


Background 


Excess UL 


Flux UL 


GeV 


cut (deg) 




events 


(95% CL) 


(cm 2 s 


-'TeV- 1 ) 


140 - 200 


8 


487 


517±23 


37 


7.2 > 


IO" 11 


200 - 280 


1(1 


683 


681 ± 27 


95 


5.1 > 


10~ n 


280 - 400 


8 


254 


242 ± 16 


75 


2.2 > 


10~ n 


400 - 560 


6 


62 


73 ±9 


14 


2.4 > 


io- 12 


560 - 790 


4 


32 


27 ±5 


27 


2.7 > 


io- 12 


790 - 1120 


4 


4 


5.7 ±2.4 


5.8 


3.7 > 


io- 13 



4. Comparison with models 

In Figure [2] we compare our flux upper lim- 
its with the theoretical 7 -ray spectra calculated by 



Bedna rek & Si tarek ( 2007h . In this model, leptons are 
injected into the GC volume according to a power-law 
spectrum, upon acceleration in the shocks produced in 
the collisions of the pulsar winds of several MSPs. 7- 
rays are then produced via Inverse Compton scattering 
of photons from the microwave background radiation 
and the thermal radiation arising from the whole GC. 
Thus, the comparison of our experimental upper limits 
with the different theoretical gamma-ray spectra al- 
lows us to constrain the total power of injected leptons 
(L c ). For this, we require the theoretical gamma-ray 
spectra to be lower than all the obtained experimental 
limits, we take into acco unt the light field of the GC 
( Bednarek & Sitarekll2007 ) and assume a distance of 7 
pc. The upper limits to L c are reported in Table 2 for 
different assumptions of the spectrum of the injected 
leptons, i.e. for different values of the spectral index 
a between the minimum energy E m \ n and the max- 
imum energy defined by the escape of leptons from 
the shock. In the case of mono-energetic injection 
of leptons we consider two different energies, 1 TeV 
and 10 TeV. Assuming characteristic values for the 
parameters of the MSPs in globular clusters (surface 
magnetic field 10 9 G and rotational period 4 ms), we 
can translate the limits to L c into limits to the product 
of the required number of the MSPs in Ml 3 (JVmsp) 
times the efficiency of the rotational energy conver- 
sion of MSPs into relativistic leptons (77), sh own in 
Table [H For example, in the case of M13, Tavanil 



1993) predicts the existence of 100 MSPs. On the 



other hand, the efficiency of lepton injection from the 
inner magnetospheres of millisecond pulsars has been 
estimated to be 77 ~ 0.1 in terms of the ex tended polar 
gap model by Muslinov & Hardingl l 1997 ). Therefore, 
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Fig. 1 . — Distribution of Alpha for the selected 7-ray 
candidates from the source (black dots) and the back- 
ground (histogram) regions. 
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Fig. 2.— The MAGIC 7-ray flux upper limits for Ml 3 
compared with spectra expected for the range of pa- 
rameters of the model sho wn in Figure 9 and 10 of 
Bednarek & Sitarekl (120071) . The specific 7-ray spec- 
tra are calculated for lepton upper energy cut-off at 3 
TeV and lower energy cut-off at 1 GeV (black thick) 
and 100 GeV (black thin), and power-law spectral in- 
dices of 2.1 (solid) and 3 (dashed). The 7-ray spectra 
produced by mono-energetic leptons of 10 TeV and 1 
TeV are shown by a grey solid curve and a dashed one 
respectively. All calculations are computed assuming 
the conservative value of 1 for the free parameter of the 
model A^msp 1 V- The Whipple differential upper limit 
sh own here has bee n derived from the integral quoted 
in Hall et all d2003l) assuming a spectral index of 2.6. 
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the likely value of the product, Nmsp ' ?7> should be 
of the order of ~ 10. We show in the correspond- 
ing row of Table [2] our estimate of the upper limit to 
this product for different models of injected spectra 
of leptons. Moreover, our limits in this product are at 
the same level than the recently published ones by the 
HESS collaboration ma king use of data from 47 Tuc 
jAharonianetalJl2009h . Figure [3] shows these limits 



in the Nms p , V plane, such that for each set of model 
parameters the area above the corresponding curve 
is excluded at 95% CL. For most of the considered 
models TVmsp ' V is significantly below ~ 10. The 
only exception is the model with the steep spectrum 
of leptons which extends down to 1 GeV. Note that 
even if the number of M SP in Ml 3 is only equ al to 
5 (as presently observed, Camilo & Rasiol 2005b . we 
can already obtain the acceleration efficiency of lep- 
tons to be ~0. 1 in the case of their injection with the 
hard (spectral index 2. 1) and mono-energetic spectrum 
respectively (see Figure[3]). 

The 7-ray spectra produced in the curvature pro- 
cess inside the inner pulsar magne tospheres are pre- 
dicted to cut-off below ~ 100 Ge V dBulik et alJ koOO: 
Harding. Usov & Muslinovll2005l) . Thus, they can not 
extend to the energy region investigated by our mea- 
surement. An Inverse Compton 7-ray component is 
expected from leptons accelerated in the inner magne- 
tospheres which extends > 100 GeV. But its flux is 
predicted to be at the level of several orders of magni- 
tude below the here presented upper limits for most of 
the energy range c overed by our limits. For example, 
Venter el al. ( 20091) computation of VHE spectrum for 
GC Tuc 47 and Ter 5 pre dicts a flux level similar to 
the ones of the model by iBednarek & Sitarekl {2007) 
for these pulsars only for narrow energy band above 
E = 1 TeV. 

Therefore, we conclude that the inner magneto- 
sphere 7-ray emission of millisecond pulsars is not 
likely to be detected by present observations with an 
analysis threshold of the order of ~ 100 GeV even 



Table 2: Upper limits on the power of injected leptons 
(L c ) and in 7V M sp • V 



Emin 


100 GeV 


100 GeV 


1 GeV 


1 GeV 


mono: 


mono: 


a 


2.1 


3.0 


2.1 


3.0 


1 TeV 


10 TeV 


L„ 


0.6 


1.0 


1.0 


60 


0.2 


0.5 


(xlO^ergs" 1 ) 














Nmsp ■ V 


0.5 


1.0 


1.0 


50 


0.2 


0.4 



from globular clusters containing hundreds of MSPs. 

5. Conclusions 

We have obtained the strongest upper limits to date 
on the VHE 7-ray flux from the massive globular clus- 
ter Ml 3. Our upper limit is ~ 2 times lower than 
the previous limit for VHE energy emission from Ml 3 
quoted by WHIPPLE, and extends to energies down to 
140 GeV. Our upper limits allow us to constrain the 
population of the millisecond pulsars expected in Ml 3 
and the acceleration scenarios of leptons by millisec- 
ond pulsars. Our result strongly suggests that either the 
number of millisecond pulsars in Ml 3 is significantly 
lower than the estimate of ~ 100, or the energy conver- 
sion efficiency from millisecond pulsars to relativistic 
leptons is significantly below the value quoted in re- 
cent modeling of high energy processes in the magne- 
tospher es of millisecond pul s ars. O ur upper limits re- 
garding |Bejlmu£k_&_Sjtarek| (120071) model parameters 
are the same level than the ones obtained by the HESS 
collaboration making use of observations of Tuc 47. 
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Fig. 3. — Exclusion contours in the Nmsp, V plane for 
the different models considered in the text. The model 
with parameters a=3 and cut-off at 100 GeV overlaps 
the a=2.l with cut-off at 1 GeV one. The horizon- 
tal and vertical long-dashed black lines show the refer- 
ence values of 77=0. 1 and Nmsp = $ respectively. 
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